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Abstract 

In this paper I discuss the extension of the analogy between gravita- 
tion and some systems of condensed matter physics from kinematics to 
dynamics. I will focus my attention on two applications of the analogy 
to the dynamics of fluids that have been recently proposed: the study of 
backreaction effects and the calculation of the depletion in Bose-Einstein 
condensates, showing how this extension is possible and stressing the main 
differences with respect to the gravitational context. I will conclude with 
some remarks about the actual reliability of the proposed scheme, pointing 
out the basis issues that have still to be addressed. 

One of the most peculiar features of actual theoretical physics is that some of 
its branches suffer an incredible lack of experimental data that seriously risks to 
compromise their development. Among them is Quantum Field Theory (QFT) 
in curved space, with its most famous result - the Hawking effect lj, widely 
considered a milestone of theoretical physics - still unobserved. In contrast with 
what previously believed about black holes, Hawking, studying the propagation 
of quantum fields on the background of a collapsing body forming a black hold3, 
showed that a thermal flux of particles emerges after the horizon forms, and 
therefore the black hole radiates. It is of fundamental importance that it is a 
pure kinematical effect, related only to the geometrical propagation of the field 
modes next to the horizon and completely independent of the details of the 
collapse itself. Therefore it will occur whenever the same spacetime structure 
emerges. 

As already pointed out, for the extreme weakness of the produced effect, the 
Hawking radiation has so far no experimental observation confirming it: in fact 
for astrophysical black holes the temperature associated to the Hawking emis- 
sion is of order 1CP 8 K or lower, far beyond the CMB temperature (~ 3 K) 

1 This effect has been obtained within the semiclassical theory of Gravity, where the space- 
time is described by the classical Einstein theory of General Relativity (GR), while the fields 
propagating on this background are quantized. 



which hoplcssly covers it. A possible test in next generation accelerators colli- 
sions has been also proposed, provided an ad hoc extra-dimensions theory lowers 
the scale of Quantum Gravity at ~ TeV [2J. 

A more promising way out can be in the so called Analog Models, i.e. some 
systems of condensed matter which, under some circumstances, behave as grav- 
itational ones. As shown by Unruh in his seminal work [3J, the fluctuations 
in an Eulerian fluid propagate exactly as a massless scalar held on a fictitious 
curved background, whose metric tensor depends on the characteristics of the 
underlying fluid. In fact an inviscid irrotational fluid is described by the hydro- 
dynamical equations 

dtp + v • { P v) = o , d t e + x -v 2 + p{ P ) = o , (i) 

with p the density, v = V0 the velocity, p(p) = *p and u the internal energy. 
Expanding eqs.([T]l at first order in the linear fluctuations p\ and 0\ around a 
mean field solution of eqs.([T]) one gets 

□01=0, /5x = --§(^0i+V0-V0i) , (2) 

where □ = -^^d^^J—gg^dv) is the covariant dAlembertian calculated from 
the so called "acoustic metric" 

p ( -(c 2 -^) -if? \ 
9^=-[ ! J- ( 3 ) 

and c is the sound speed c 2 = p 

9 Pa 

Thus the field 0\ satisfies the equation of motion for a massless scalar field in 
the curved background described by the tensor metric ([3J. For a homogeneous 
fluid the acoustic metric reduces to the flat minkowskian one, while it can show 
a black hole form if \v\ — c: there goo vanishes and a sonic horizon forms [J]. 
The supersonic region is the acoustic analog of the black hole, where sound is 
trapped as light is in the gravitational field of a black hole. 
Once the hydrodynamical fluctuations are quantized, following Hawking's ar- 
guments, a flux of phonons (the sonic quanta) at a temperature Th = g^ntsc 
is expected, where k — ^V(c 2 — v 2 )\h or is the acoustic analog of the surface 
gravity of the horizon [3J. This is the sonic analog of the Hawking effect. For 
this reason hydrodynamical systems can in principle be used as "laboratories" 
to observe this effect (among other effects otherwise unfathomable) . One of the 
most promising setting to this aim are Bose Einstein condensates (BEC), for 
which Th ~ 10 nK with respect to a condensate temperature T c w fOO nK 0. 
Unfortunately, the analog Hawking radiation is still difficult to point out di- 
rectly, both for the fact that Th is actually very low and for the noise produced 
by other effects that risks to cover it. Yet the recent technological improvement 
in handling BEC and the possibility to study analog models more in depth make 
us confident that it can be soon observed. 
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Dynamics 



As shown in the previous section, the analogy works perfectly at the kinematical 
level: it occurs for the amazing fact that sound waves in a fluid "feel" inho- 
mogeneities exactly as a quantum field on a gravitational background "feels" 
curvature and it is therefore possible to use the GR formalism to describe their 
motion. But the features of the acoustic metric strictly depend on the under- 
lying medium and its flow, which can be described in no way using the GR 
framework: they are indeed non-relativistic (c SO md <C cu g ht) and newtonian. 
The equations governing the dynamics in the two setting evidently differ: on 
one side are Einstein equations, on the other hydrodynamical equations. 
Does this mean a full stop in the application of analogy beyond the kinematical 
issues? Or is there still the possibility to go further addressing also dynamical 
ones? 

I will show how is still possible, with special care, to use analog models to show 
up also behaviors depending on the dynamics. Of course the analogy will be not 
"perfect" any more and the outcomes in the condensed matter context could 
have nothing to do with what known in GR. Yet the formalism developed in 
QFT in curved space turns out to be a powerful tool to evaluate analytically 
quantities otherwise nasty to evaluate with the standard techniques. 
As examples, I will show in the following two seemingly very different applica- 
tions of this formalism: the backreaction of the analog Hawking radiation and 
the depletion in BEC. 

Backreaction in sonic black holes 

As already mentioned, Hawking radiation is still difficult to be directly observed 
even in the acoustic setting. It could be easier to try and see it indirectly 
through the modifications its presence will induce on the fluid: this beck-effect 
is the so called backreaction. Of course it will depend on dynamics and there is 
no reason to expect that what found for gravitational black holes still holds for 
these completely different objects. 

In GR backreaction is governed by the semiclassical Einstein equations, namely 
Gp,v{gafi) = Ji(Tfj, v (ga/3)}, where (T^) is the stress energy tensor of the quantum 
field which acts as a source in the semiclassical equations and modifies the 
background spacetime. The first issue is to find the analog of the semiclassical 
equations for the hydrodynamical systems. In [5] it is shown how to solve this 
problem going further in the expansions of eqs.([T]) up to quadratic order in the 
fluctuations, yielding 

dtpB+y-{pBV B ) = -di(V^{T 0i )) , i]jb + Iv1+Kpb) = \^-{T) , (4) 

I I c 

with the subscript b meaning " corrected by backreaction" . 
In eqs.(UJ) the pseudo-stress energy tensor for the fluctuations (namely for a 
massless scalar field 6\) has been introduced: 7) t „ = d^Oid^Ox — \g^ v d a Q\d a Q\ 
and T is its trace. This procedure is equivalent to what one usually do in QFT in 
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curved space to write down the backreaction equations for gravitational waves. 
It is worth noting that eqs.([4]) are the most general equations governing backre- 
action effects in the hydrodynamical setting for c = const: they are not linked 
to the presence of a sonic black hole and describe classical as well as quantum 
fluctuations backreaction. All the informations relative to the different situa- 
tions are enclosed in (T^ v ): the brackets are stochastic mean values in describing 
classical fluctuations or expectation values in the quantum case. The quantum 
state which properly describes black hole formations and Hawking radiation is 
the so called Unruh state [7]. 

Under some simplifying hypothesis, eqs.(g]) has been solved analytically [B] and 
numerically [8| for a sonic black hole in the Unruh state. It has been shown 
that, due to the presence of the Hawking radiation, next to the horizon the fluid 
slows down and the black hole size shrinks. Moreover the emission temperature 
gets lower and lower as the evaporation goes on, reaching zero temperature 
emission in an infinite amount of time. This is not a surprising behavior: the 
evaporation proceeds at the expenses of the kinetic energy of the fluid, which 
in fact slows down, and the temperature (related to the velocity on the hori- 
zon) also decreases. Yet it is completely different from what is known about 
Schwarzschild black holes, that get hotter and hotter as the evaporation goes 
on. The sonic black hole looks rather like a near-extremal charged black hole: 
this kind of black hole reaches zero temperature emission in an infinite time, 
but with different power law with respect to the sonic ones [5]. 

Depletion in Bose-Einstein condensates 

As well as other ultra-cold atoms systems, BEC can be described in some regimes 
by hydrodynamical equations [lOj . In the mean field limit the bosonic operator 
^ describing the system can be split as ^ = ip + ip: ip = describes the 
condensed part, while <p is a quantum fluctuation. Under some hypothesis 
and writing ip — ^fpe l6 / h , the condensed part satisfies the classical equation of 
motion analog to eqs.fl]), where now p is the density of particles and v = m _1 V0 
the velocity of the fluid. In practice this description is possible if the scale 
of variation of the inhomogeneities L is much bigger than the healing length 
(, = %/ (V2mc), that is related to the mean distance between the constituents of 
the fluid. Moving to the fluctuation part, without loss of generality, ip can be 
rewritten in the acoustic representation as 



with p\ and 9\ two real fields. The equations describing them assume the 
hydrodynamical form and can therefore be rewritten using the gravitational 
formalism as in eqs.© [4]. It is worth stressing that this description is correct 
only for modes with wavelength A > (: the modes with very high frequency 
in fact do not obey the relativistic dispersion relation lo = ck implicit in the 
GR framework, but the full Bogoliubov dispersion relation for BEC, namely 
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U) = ckyjl + ^k 2 . 

In the hydrodynamical regime it is therefore possible to use BEC as labs for test- 
ing GR effects, but also to handle the analogy to give analytical estimations of 
quantities that are hardly evaluable with the standard tools. In this spirit, QFT 
in curved space has been used to work out the depletion for a non-homogeneous 
BEC [TT] . The total particles density can be written as a sum of condensed and 
non condensed part 

Ptot = (* t *) =ip*ip+($<p) =p + p. (6) 

The last term describes the so called depletion: because of quantum and tem- 
perature fluctuations, a non vanishing fraction of the total particles number is 
not in the condensed phase. The depletion is explicitely known only for homo- 
geneous BEC where p = p = const. At low temperature (kbT <C Pog), it reads 

\M 



P° = ^Tj2 (A)<*) 3/2 



V 2 Po9 ) 



(7) 



Moving to non homogeneous condensates with vanishing velocity, the depletion 
is usually approximated by eq. ([7]) with now po replaced by the actual space-time 
varying density p(x, t) , in the so called local density approximation (LDA) [12] . 
It is clear that this approximation is not able to take the effects of spatial deriva- 
tives into account. It is worth noting that the derivation of the zero temperature 
quantum depletion, which does not vanish in the T — > limit, is based on the 
full Bogoliubov dispersion relation for the modes. Unlike this, the low tempera- 
ture correction comes from the phononic part of the spectrum only, i.e. that one 
assumed in the GR analogy. The integral giving the thermal depletion is domi- 
nated by modes with wavelength A 3> £, so the result is insensitive to the high 
k behavior of the dispersion relation. Therefore, while the correct description 
of the ground state requires the knowledge of the full "microscopic" theorjH, 
deviations from the ground state value can be evaluated using the low energy 
"macroscopic" phonon theory. In other words, we assume that the corrections 
coming from the inhomogeneities of the condensate to the ground state result, 
sensitive to the whole quantum theory describing the system, are determined by 
those modes with wavelength comparable with the scale of the inhomogeneities 
(as curvature is seen by modes with A of the same order of the curvature it- 
self). Those modes therefore satisfy the hydrodynamical approximation request 
(A ~ L ^> £) and for them the analogy works. □ 

Using techniques developed within QFT in curved space, it is possible to over- 
take the difficulties in handling analytical calculations for non-homogeneous 
BEC pj]. Now, noting that 

Pm = {LpLp) = ^o- + — *-fT\ 1 + °{t 2 )) ' (8) 



2 The linear phonon dispersion relation, characteristic of a massless relativistic free field 
would lead (after renormalization) to a vanishing quantum depletion. 

3 This is more than a guess, since it is partially confermed by the results in this direction 
coming from the Casimir effect [13] ■ even if so far there is not a conclusive answer. 
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the calculation of the depletion has been reduced to the calculation of the ex- 
pectation value of the squared of a massless scalar field minimally coupled to the 
acoustic spacetime ([3]). For a spherically symmetric static BEC configuration, 
using renormalization techniques peculiar of QFT in curved space, it is possible 
to give an explicit description for the thermal and quantum depletion for an 
inhomogeneous BEC in term of a general density p(r): 
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(9) 



where m is the mass of the atoms and g is a constant proportional to the scat- 
tering length a: g = 4irh 2 a/m. Note that the thermal term exactly reproduces 
the thermal depletion in LDA. The other terms instead generalize the zero- 
temperature depletion in LDA taking the spatial derivatives of the condensate 
density explicitely into account. In fact all of them vanish for a homogeneous 
BEC. 

Moving to configurations with a radial stationary velocity (v = v(r)r/\r\), the 
concept of thermal equilibrium loses its meaning, but it is still possible to eval- 
uate the effects of this velocity on the quantum depletion, yielding [TT] 



96ir 2 hc 2 
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(3c 2 c' - vV - 2czw')(3c 2 c' - 2v 2 c' - cvv') 



(10) 



where we have introduced the local sound speed c, that for BEC is related to the 
density by the equation of state mc 2 = gn. This expression generalizes eq.([9]) 
considering also contributions due to the velocity of the condensate. It reduces 
to eq.Q only for v(r) = 0: for v = const in fact a curvature is induced since a 
radial velocity can not be trivially eliminated by a Galileo transformation. 
In the case of an acoustic black hole, the expression in ea. pUl) diverges on the 
acoustic horizon (v 2 = c 2 ): this is a familiar divergence in GR, related to the 
fact that the quantum state over which the expectation value of (0 2 ) has been 
calculated is not the right one to describe the region next to the horizon. In GR 
it has been shown that for a Schwarzschild black hole the same calculation for 
the Unruh state leads to a regular result: the same is expected also for acoustic 
black holes. 



6 



Conclusions 



In this paper I have illustrated the first step in extending the analogy between 
gravitation and condensed matter systems also to dynamical issues. It has been 
stressed how this extension has to be handled with particular care, without 
expecting in the acoustic setting an exact counterpart of the results already 
obtained in GR (unlike what happens with the kinematical ones). 
However there is still a crucial issue to be addressed: how big are the errors in 
neglecting high momenta effects? 

This question involves both QFT in curved space and condensed matter physics. 
As already mentioned, for BEC as in general for fluids, very high k modes see 
the coarse grain of the medium, and, as a consequence, they do not obey the 
relativistic dispersion relation, but a short-distance modified one. 
In the gravitational context instead the question is related to the so called 
transplanckian problem: the frequency of the modes which give the Hawking 
radiation is indeed transplanckian, i.e. beyond the actual reliability limit of 
the theory used to study the process, namely QFT in curved space. This fact 
does raise doubts on the reliability of the Hawking effect itself, as a consequence 
of the breaking of QFT in curved space beyond Planck energy scale, where 
Quantum Gravity effects should become important. Thought a full Quantum 
Gravity theory is not known and therefore it is not possible to safely say how 
these behaviors will emerge, it is widely believed that a first hint of such effects 
should show up in a short-distance modification of the dispersion relation [14) 
(like what happens in the acoustic setting). Recent investigations on Hawking 
effect within a QFT with modified dispersion relation have shown that Hawk- 
ing radiation is only slightly affected by these modifications, turning out to be 
unaltered in its fundamental features [T51 [§] • 

This is a good starting point to carry out deeper investigations involving also 
the calculation of expectation values in presence of such modified dispersion 
relations. Moreover also the first result concerning the calculation of renormal- 
ized expectation value goes in this direction: in fact the Casimir energy density 
between two conducting plates can be evaluated with great accuracy using the 
phonon dispersion relation, provided the plate separation is much larger than 
the interatomic distance of the constituents of the plates [13) . Similar analysis 
has been performed also in the cosmological context, yielding so far only partial 
results [16] . Nothing has still been obtained for a black hole background. 
A deeper knowledge of these effects will be of crucial importance in several field: 
cosmology, black hole physics and cold-atoms systems. A more accurate use of 
QFT in curved space will turn out to be very useful both in studying gravitation 
directly and in leading a more reliable use of the analogy in the dynamics of the 
condensed matter systems. 
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